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We investigate the potential for the Deep Underground Neutrino Experiment (DUNE) to probe 
the existence and effects of a fourth neutrino mass-eigenstate. We study the mixing of the fourth 
mass-eigenstate with the three active neutrinos of the Standard Model, including the effects of new 
sources of CP-invariance violation, for a wide range of new mass-squared differences, from lower 
than 10“® eV^ to higher than 1 eV^. DUNE is sensitive to previously unexplored regions of the 
mixing angle - mass-squared difference parameter space. If there is a fourth neutrino, in some 
regions of the parameter space, DUNE is able to measure the new oscillation parameters (some very 
precisely) and clearly identify two independent sources of CP-invariance violation. Finally, we use 
the hypothesis that there are four neutrino mass-eigenstates in order to ascertain how well DUNE 
can test the limits of the three-massive-neutrinos paradigm. In this way, we briefly explore whether 
light sterile neutrinos can serve as proxies for other, in principle unknown, phenomena that might 
manifest themselves in long-baseline neutrino oscillation experiments. 

PACS numbers: 14.60.Pq, 14.60.St 


I. INTRODUCTION 

It is now established, beyond reasonable doubt, that neutrinos have mass and that leptons mix. In order to further 
explore the neutrino sector and exploit the oscillation phenomenon, ambitious next-generation, long-baseline neu¬ 
trino oscillation experiments are under serious consideration, including the Deep Underground Neutrino Experiment 
(DUNE) proposal in the United States [1] and the HyperKamiokande detector (and accompanying neutrino source in 
J-PARC) in Japan [5]. The goals of these projects include the search for leptonic CP-invariance violation and testing 
the limits of the so-called three-massive-neutrinos paradigm (for more see, for example. Ref. [5]). 

The hypothesis that there are three neutrinos, at least two of them massive, and that these interact as prescribed 
by the Standard Model of electroweak interactions, accommodates almost all neutrino data. There is, however, plenty 
of room for new phenomena. The unitarity of the leptonic mixing matrix, for example, has not been thoroughly 
explored* (see, for example. Ref. S)- In spite of tremendous experimental progress, little nontrivial information 
regarding the appropriateness of the three-massive-neutrinos paradigm has been collected over the past seventeen 
years. 

Next-generation long-baseline experiments can probe several different new phenomena, including new, weaker-than- 
weak, interactions involving neutrinos and charged-fermions [S] that lead to anomalous matter effects, new long-range 
forces [6], and the existence of very light new states [7j. Here, we concentrate on the simple hypothesis that there is 
a fourth neutrino mass-eigenstate of unknown, but very small, mass - 7714 less than a few eV - and assume that there 
is a nonzero probability that the new neutrino state can be measured as one of the active neutrino states, i.e., we 
assume the leptonic mixing matrix U to be 4 x 4 and that Uai, a = e, ij,,t are nonzero. Electroweak precision data 
require the fourth neutrino flavor-eigenstate not to interact with the W and Z bosons with Standard Model strength, 
so we refer to it as a sterile neutrino. 

We concentrate on this hypothesis for a few different reasons. It is simple and easy to parameterize, and very 
familiar [7]- Indeed, certain aspects of the effects of sterile neutrinos on different long-baseline experiments have 
been studied in the recent past (see, for example. Refs. [SHIl]). Sterile neutrinos are also a very natural and benign 
extension of the Standard Model and could be a side effect of the mechanism responsible for the nonzero neutrino 
masses (see, for example. Refs, [inilllj). Furthermore, there are the so-called short-baseline anomalies [T7IET| . These 
might be pointing to more new physics in the leptonic sector, but a convincing, robust explanation remains elusive. 
Sterile neutrino interpretations to the short-baseline anomalies are, arguably, the simplest explanations of these data. 
Here, we remain agnostic regarding the new-physics origin of the short-baseline anomalies but, on occasion, will 
highlight the region of mass and mixing space that is preferred by them. Finally, we would like to explore whether 


* Arguably, the data only allow unitarity checks for the first row, 2 3 IHeiP 


1, and third column, Ua=e,M,T IHoSp = 1- 
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light sterile neutrinos can serve as proxies for other, in principle unknown, phenomena that might manifest themselves 
in long-baseline neutrino oscillation experiments. This aspect of our analysis will become more clear later. 

This paper is organized as follows. In Section|nj we review four-flavor neutrino oscillations. Since we are interested 
in a large range of new mass-squared differences, we pay special attention to the neutrino oscillation probabilities in 
the limits when the new oscillation frequency is much larger or much smaller than the known oscillation frequencies. 
In Section III we discuss in detail the capabilities of DUNE to (a) rule out the sterile neutrino hypothesis assuming the 
data are consistent with the three-massive-neutrinos paradigm; (b) determine the new mixing parameters assuming 
there is one sterile neutrino. Here, we assume new mass-squared differences that range from 10“® eV^ to 1 eV^; and (c) 
diagnose that there is physics beyond the three-massive-neutrinos paradigm assuming there is a fourth mass-eigenstate. 
In Section [IV[ we offer some concluding remarks. 


II. OSCILLATIONS WITH FOUR NEUTRINOS 


We consider a fourth neutrino Vg that does not participate in the weak interactions but that can mix with the 
other three neutrinos of the Standard Model. The misalignment between the mass-eigenstates Vi {i = 1, 2, 3, 4) and 
flavor-eigenstates Va (ex = e, /i, r, s) can be described by a general 4x4 unitary matrix parameterized by six angles 
(j)ij {i, j = 1, 2, 3, 4; i < j) and three phases iji, r]2, rjs. We choose the matrix elements Uai to be:^ 


Ue2 = S 12 C 13 C 14 , 

Ue 3 = e“*''^Si3Ci4, 

17e4 = e“*''^Si4, 

Uf_i2 = C24 (ci2C23 — e*''^Si2Si3S23) ~ ’^^^Sl2Sl4S24Cl3, 

Ufj.3 = S23C13C24 ~ ^'^^Si3Si4S24) 

Ufj,4 = e“*''^S24Cl4, 

Ut 2 = C34 ( —C12S23 — e*''^Si2Si3C23) — 6*’'^ C13C24S12S14S34 
— (ci2C23 — e*’'^Si2Sl3S23) S24S34, 

Ut3 = C13C23C34 — ’^^^Si3Si4S34C24 — e*’’^S23S24S34Cl3, 

Ut4 = S34C14C24, 


(11.1) 

(11.2) 

(11.3) 

(11.4) 

(11.5) 

( 11 . 6 ) 

(11.7) 

(11.8) 

(11.9) 


where Sij = sin (j)ij , Cij = cos (j)ij ■ The matrix elements not listed here can be determined through the unitarity 
conditions of U. 

When the new mixing angles ^ 14 , (f) 24 , and 4)34 vanish, one encounters oscillations among only three neutrinos, and 
we can map the remaining parameters { 4 >i 2 , 4 'i 3 i 4 ’ 23 i Vi} {^ 12 ) ^i 3 j ^ 23 j ^cp}-: the well-known mixing parameters 
that define the standard 3x3 leptonic mixing matrix in the three-massive-neutrinos paradigm (using the Particle 
Data Group convention [22]). In the limit where (j)i 4 , (j> 24 , and 4>34 are small, the angles (j)i 2 , 4>i3 and (j )23 play roles 
very similar to those of 0 i 2 , 9 13 and 023! respectively. We discuss this in more detail in Section |ITT| The best-fit 
values from Ref. [12] for a three-flavor fit to existing data are sin^ 0 i 2 = 0.308 ± 0.017, sin^ 0 i 3 = 0.0234tQ gQ^g, and 
sin^ 023 = 0.437lg;g23; the CP-odd phase 6 cp is virtually unconstrained. 

The amplitude for Va to be detected as Vjs after propagating a distance L in vacuum is 


Aap = + U^2U*P2 - 1) + Pa3C^3 - 1) + C7a4C^4 - 1), 


(II.IO) 



^ We ignore potential Majorana phases because they do not affect oscillations in any realistically observable way. 

^ If 1^4 decoheres from the other three neutrinos, then the expression for is modified by neglecting the interference of the oscillations 
related to A14 with those of A12 and A13, cf. Eq. ( |II. 12 [ |. Decoherence occurs if, for example, U4 is produced incoherently, or, during 
propagation, the U4 wavepacket becomes well-separated from the wavepacket containing U2 ^-nd 
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eV^, Ami 3 = 2.43 x 10 ^ eV^, assuming the neutrino 
The value of 77 x 4 is mostly unconstrained, so can 

i.e., 


the mass-squared splittings Am 42 and Am\^ are very close to the ones that fit the neutrino data assuming there are 
only three neutrino species (explicitly, Am \2 = 7.54 x 1 0~^ 
mass hierarchy is normal [H]), as we will discuss in Sec. Ill 
be larger or smaller than Amf 2 and Am\^. We do, however, restrict our analyses to positive Amf^, i.e., 1714 > mi. 
In summary, including the fact that we will always assume the normal neutrino mass hierarchy for the mostly active 
states, our masses are ordered as follows: toi < m 2 < m 3 , and 7714 > mi. As we vary Am^^ from very small to very 
large, we allow for all different mass orderings: 7774 < m 2 < m 3 ; m 2 < m 4 < m 3 ; and m 2 < m 3 < m 4 . 

The amplitude simplifies considerably when A 14 <C 1. In this limit, the last term in Eq. (II.10) is small compared 
to the others, so 


Pal3 — 


Scp + U^2U;2 - 1 ) + U^3U;3 - 1 ) 


( 11 . 11 ) 


Because the experimental normalization uncertainties we will consider are 0(1%), the oscillations associated with 
Amf^ will not be discernible if A 14 < 10“^ over the entire range of reconstructed neutrino energies. For long-baseline 
oscillation experiments with L ^ 0(10^ km) and E^, ~ 0(1 — 10 GeV), this condition translates into Am \4 < 10“^ 
eV^.§ Nonetheless, in this scenario, oscillations can be distinct from those among only three neutrinos. Here, the 
elements U^i, a = e, t; 7 = 1, 2, 3 do not form a unitary matrix and the number of independent parameters is larger 
than four, including sources of OP-invariance violation beyond the phase rji [nin]. We return to this in Sec. 

When Ai 4 ^ 1 but z /4 is light enough to be produced coherently in the initial neutrino state, the oscillations of 
Am 34 can be too rapid to be resolved by the finite energy resolution employed by the experiment. The oscillations 
of Am 34 average out if A 14 x (SE/E,^) is, roughly, larger than 27r, where the energy bin width is SE and the bin’s 
central energy is Pjy. For L ~ 0(10^ km), ~ 0(1 — 10 GeV), and SE ~ 0.25 GeV, this will occur if Am 44 > 1 
eV^. In this case, 


III 


PaP — 


Scp - + U^2 U*p2 - 1) + " l) 


Pa4P, 


/34| 


( 11 . 12 ) 


This limit is, as far as measurements of the oscillation probabilities are concerned, equivalent to the decoherence of 
1/4 from the other neutrinos. As in the A 14 <C 1 limit, oscillations are in general distinct from those among only three 
neutrinos. For example, one is, in principle, also sensitive to sources of CP-invariance violation beyond the phase rji 

[Hill]. 

When neutrinos propagate through matter, elastic, coherent, forward scattering modifies the oscillation probabilities 
in a well-known way. This can be parameterized via an effective potential generated by the background of electrons, 
protons and neutrons. The Hamiltonian SElap that describes neutrino oscillations, in the flavor basis, is |23j 


/ SHgp \ 
V1 km“^ / 



Sge.S. 


'Pe 



SgsSps: 


(11.13) 


where (A/1 eV^) = (3.85 x 10“‘*)ye(p/l g cm“^) characterizes the charged-current interactions, (A'/l eV^) = (1.92 x 
10“’^)(1 — Ye){p/1 g cm“^) characterizes the neutral-current interactions. Ye is the electron fraction for the matter 
background, and p is the density of the background.^ The signs of A and A' are flipped for antineutrinos. The 
Earth’s crust typically has 1/ ~ 0.5 and p ~ 3 g cm“^ [53]. In the presence of matter, the Hamiltonian is no longer 
diagonal in the mass basis and the exact expressions for the oscillation probabilities are much more cumbersome. In 
our analyses, we treat the flavor evolution of the neutrino states numerically. 

Our analysis makes use of initially muon-type neutrinos produced in pion decay to study P^^ and P^e- Because of 
experimental challenges involved in working with r leptons, we do not consider oscillations into Vt.§ ** Gonsequently, 


§ Note that A.rn ^2 close to this limit, i.e., the wavelengths of its associated oscillations are too long to significantly impact oscillations at 
such an experiment. Nonetheless, sensitivity to the oscillations associated with Arrijj comes from the interference with the oscillations 
due to Am^g. As we discuss in Appendix A, long-baseline experiments rely on information regarding Am ^2 ^12 from other sources, 
including solar and reactor neutrinos, in order to precise ly me asure all oscillation parameters. 

^ The zero-point of the potential has been shifted in Eq. (11.13 1 so that the neutral-current contribution appears with the opposite sign 
in the sterile—sterile part of the Hamiltonian. 

** The study of tau appearance requires neutrino energies above the tau-production threshold for neutrino—nucleon scattering, around 
3.4 GeV. Hence, for the energies under consideration here, tau-appearence is severely phase-space suppressed. Furthermore, detectors 
must be able to identify taus with nonzero efficiency, an issue that is actively under investigation. 
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we do not expect to learn much about (/) 34 , which only appears in the matrix elements Uri- For long-baseline oscillation 
experiments with L ^ 0(10^ km) and ^ 0(1 — 10 GeV), is mostly sensitivity to (^ 24 , while Pp,e is sensitive to 
both ^24 and ^ 14 , mostly via the product sin ^24 sin (/)i 4 . Therefore, we expect DUNE to have greater sensitivity to 
024 than to 014. Furthermore, these two channels both depend on the CP-odd phase 771 , as well as the combination 

Vs =V2-m- (11.14) 

In order to distinguish the effects of 772 from 773, one requires information regarding Uri, which, as just argued above, 
is unavailable in the absence of searches for t appearance or disappearance. 


III. EXPERIMENTAL SENSITIVITY TO A FOURTH NEUTRINO AT DUNE 

We investigate the sensitivity of the proposed Deep Underground Neutrino Experiment (DUNE) [T] to a fourth 
neutrino. We consider that DUNE consists of a 34 kiloton liquid argon detector and utilizes 1.2 MW neutrino 
and antineutrino beams originating 1300 km upstream at Fermilab, consistent with the proposal in Ref. [1]. We 
also simulate that the detector has resolution of a [GeV]= 0.15/i/Py [GeV] for identifying electrons and cr [GeV] = 
0.20El, [GeV] for identifying muons. The neutrino energy ranges between 0.5 and 20 GeV and the flux is largest 
around 3.0 GeV. In the following analyses, we simulate six years of data collection: 3 years each with the neutrino 
and antineutrino beams. 

We use the neutrino fluxes and signal reconstruction efficiencies projected by DUNE (Fig. 3.18 and Table 4.2 in 
Ref. [1], respectively) and the neutrino-nucleon cross-sections reported in Ref. [24] to calculate expected yields. For a 
three-neutrino scenario, we use input values consistent with the best-fit results compiled in Ref. [52]: sin^ 6*12 = 0.308, 
sin^ 6*13 = 0.0235, sin^ 6*23 = 0.437, A 77712 = 7.54 x 10“® eV^, Aml^ = -1-2.43x 10“^ eV^ (hence a normal neutrino mass 
hierarchy), and Sqp = 0 . The four dominant backgrounds are consequences of muon-type neutrino neutral-current 
scattering (‘V^ NC”), tau-type neutrino charged-current scattering (“ 77 ^ —^ Pr GC”), muon-type neutrino charged- 
current scattering (‘V^ —>■ GG”), and beam electron-type neutrino charged-current scattering (‘Ve —>■ Ve beam 
CG”), depicted in Figs. [^a)-(d). The rates associated with these backgrounds are taken from Ref. [Tj. We reproduce 
the signal and background yields in Ref. [T] for the appearance (P^e) and disappearance (P^n) channels, shown as 
dashed lines in Figs, [^a)-(d), i.e., “ 77 ^ —)■ Vg signal 37/” and “ 7 /^ —)■ 7 /^ signal 37/,” respectively. In Appendix A, we 
demonstrate comparaHe sensitivity to those computed in Ref. [Tj. 

In order to illustrate the effects of a fourth neutrino, the expected yields along with the three-neutrino yields in 
Fig. are depicted for sin^ 0 i 4 = 0.023, sin^ 024 = 0.030, sin^ 034 = 0, A 77714 = 10“^ eV^, and rjs = 0 (“ 7 /^ —>■ T/g 
signal 47 /” and “ 7 /^ —>• 7 /^ signal 47/”). We choose the value of A 77714 such that several oscillations due to the fourth 
neutrino occur within the energy window of the experiment. Here, the input values of 0 i 2 , 0 i 3 , and 023 are slightly 
different from the values mentioned above for 6 * 12 , di 3 , and 6 * 23 , and are chosen so that the values of |Ue 2 p, |Ue 3 p, 
and |U;i 3 p are consistent with three-flavor fits to the neutrino data [ 22 ].* 


A. Constraining the Four-Neutrino Hypothesis 

If the data are consistent with the three-neutrino hypothesis - the three-neutrino scenario outlined above - one can 
place upper bounds on the values of 0 i 4 and 024 for given values of A 777 34 . We calculate 95% confidence level (GL) 
exclusion limits for a fourth neutrino in the sin^ 0 i 4 - A 77 i ^4 and sin^ 024 - A 777 f 4 planes, depicted in Fig. using the 
appearance and disappearance channels and assuming running for three years each with the neutrino and antineutrino 
beams. We include normalization uncertainties of 1% and 5% for the signal and background yields, respectively. 

Fig. j^a) also depicts the results from the Daya Bay [53] and Bugey [5^ experiments in the sin^ 0 i 4 - Aml,^ plane. 
The existing experiments have greater sensitivity to sin^ 0 i 4 for values of A 77734 > 10“"^ eV^. For smaller values of 
A 777 34 < 10“^ eV^, none of the experimental probes, including DUNE, can “see” the very long new oscillation length. 
Nonetheless, since DUNE measures both appearance and disappearance, we can constrain very large values of sin^ 0 i 4 
as these render the upper-left 3x3 mixing sub-matrix unacceptably non-unitary. The same phenomenon can be 
observed in the sin^ 024 - Am\,^ plane (Fig. [^b)). Very large values of sin^ 024 are ruled out, even for very small 
values of A 77744 . 


* We do not have the freedom to set all nine matrix elements Uai = e, /r, r; 2 = 1, 2, 3) equal to their three-neutrino best-fit values. 
Explicitly, we choose |f/e 2 p = 0.301, |f/e 3 p = 0.023, and = 0.427. 
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FIG. 1: Expected signal and background yields for six years (3y u + ‘ij V) of data collection at DUNE, using fluxes projected 
by Ref. [1], for a 34 kiloton detector, and a 1.2 MW beam, (a) and (b) show appearance channel yields for neutrino and 
antineutrino beams, respectively, while (c) and (d) show disappearance channel yields. The 'iv signal corresponds to the 
standard three-neutrino hypothesis, where sin^ 0i2 = 0.308, sin^ 0i3 = 0.0235, sin^ 023 = 0.437, Ami 2 = 7.54 x 10“® eV^, 
Ami 3 = 2.43 X 10“^ eV^, &cp = 0, while the Av signal corresponds to 4‘\2 = 0.315, sin^ i/iis = 0.024, sin^ 023 = 0.456, 
sin^ 014 = 0.023, sin^ 024 = 0.030, Amf 4 = 10“^ eV^, rji = 0, and rjs = 0. Statistical uncertainties are shown as vertical bars 
in each bin. Backgrounds are defined in the text and are assumed to be identical for the three- and four-neutrino scenarios: 
any discrepancy is negligible after accounting for a 5% normalization uncertainty. 


In the sin^ 024 - Ami 4 plane (Fig. [^b)), we also show results from the MINOS m experiment and note that 
DUNE will be sensitive to lower values of the mixing angle and the mass-squared difference, due to DUNE having 
greater expected yield and a broader range of LjE^ values. Because the disappearance channels depend strongly 
on |C/n 4 p, and have higher yields than the appearance channels, DUNE has greater sensitivity to 024 than 0 i 4 . We 
also note that, as expected and discussed in the previous section, if the mass-squared difference is either very small, 
Am .44 < 10 “'^ eV^, or very large, Ami 4 > 1 eV^, the limits are independent of the new mass-squared difference. 
These ranges correspond, respectively, to A 14 < 10“^, where oscillations due to the fourth neutrino are undetectable 
over the energy range of the experiment, and to the oscillations associated to the new frequency averaging out over 
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sin^ 014 sin^ 024 

(a) (b) 

FIG. 2: Expected exclusion limits for sin^ 0i4 (a) and sin^ 024 (b) vs. Ami 4 at DUNE (blue), assuming a 34 kiloton detector 
and a 1.2 MW beam with six years (3y u + 3y u) of data collection. The exclusion limits become independent of Ami 4 when 
the mass-squared difference is large (> 10~^ eV^) or small (< eV^). Results from the Daya Bay (red) and Bugey |26| 

(orange) are shown in (a), and results from MINOS [27] (maroon) are shown in (b). The three sets of four-neutrino parameters 
we consider in Section [III B[ listed in Table [I] are denoted by black stars above. 


the width of the energy bin at DUNE (here AE,^, the width of a bin, is 0.25 GeV). Finally, the sensitivity in the 
sin^ 024 - Atoi 4 plane extends to lower values of Arrii^ than that in the sin^ 024 - Ami 4 plane due to the higher yield 
of the disappearance channel. 

Additionally, we calculate exclusion limits at 95% CL in the 4 |t/e 4 p|C/ 4 , 4 p - Amf 4 plane^ in order to compare the 
DUNE sensitivity to the proposed short-baseline experiment z/STORM |32| , and the current long-baseline experiments 
MINOS, OPERA, and ICARUS |27l|30l|3T|. Additionally, we include the results of the global fit to all neutrino data 
assuming a four-neutrino scenario, reported in Ref. [33] . This fit includes data from short-baseline experiments, 
including the short-baseline anomalies discussed earlier. In Fig. we see that DUNE is sensitive to lower values of 
Am 44 than any existing or proposed experiment due to its access to a wider range of LjEi, values. 


B. Measuring the New Mixing Parameters 

Assuming the existence of a fourth neutrino, we explore the capability of DUNE to measure the new mixing angles 
and mass-squared difference. We choose three sets of parameters, listed in Table and denoted by black stars in 
Figs. and and calculate expected yields for the appearance and disappearance channels. Case 1 is consistent with 
the global four-neutrino fit performed in Ref. [33] (the red ellipse in Fig. [^. Here, Am\^ is large enough that we 


^ In the limit when oscillations due to the new mass-squared difference are dominant, the appearance channel oscillation probability takes 
the simple form 

- 4|t/e4p|UM4p sin^ ^ sin^ ■ (^-1) 

The effective mixing angle 9e^ is commonly used in the literature for Ue short-baseline appearance searches (see, for example, 

[ 23011 ). 
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10 “® 10 “^ 10 “^ 10 “^ 10 “^ 




Ii4\ 


FIG. 3: Exclusion limits in the 4|I/e4|^|?7^4|'* - Ami 4 plane for various existing and proposed neutrino experiments. Expected 
exclusion limits are shown at 95% CL for the proposed DUNE (blue) and i/STORM |31] (purple) experiments. Results are 
shown at 90% CL for the MINOS and Bugey [27] (orange), OPERA |30| (teal), and the ICARUS [31] (dark blue) experiments. 
Additionally, the fit to the 3 + 1 scenario including the short-baseline anomalies, reported in Ref. [33], is shown. The three sets 
of four-neutrino parameters we consider in Section [III B[ listed in Table [Ij are denoted by black stars above. 


expect the oscillations associated to the new (very short) oscillation length to average out at DUNE. Case 2 uses the 
same mixing angles as Case 1, but with a lower value of Amf 4 . The parameters are within the the reach of DUNE, but 
outside the reach of current and proposed short-baseline experiments.^ Here, Arrii^ is small enough that we expect 
the oscillations associated to the new oscillation length to be visible at DUNE. Case 3 has a much lower mass-squared 
difference, Arrii^ = 10“® eV^, but has large values of (j)i 4 and (j) 24 -^ Here, Aml^ is too small to be seen at DUNE. 
Nonetheless, as discussed earlier, the new mixing angles are large enough that nontrivial information on and ^24 
can be extracted. 



sin^ 014 

sin^ 024 

Am\^ (eV^) 

Vs 

sin^ 012 

sin^ 013 

sin^ 023 

Am\2 (eV2) 

Amf^ (eV^) 

Vi 

Case 1 

0.023 

0.030 

0.93 

—7r/4 

0.315 

0.0238 

0.456 

7.54 X 10-5 

2.43 X 10-5 

7r/3 

Case 2 

0.023 

0.030 

1.0 X 10-2 

—7r/4 

0.315 

0.0238 

0.456 

7.54 X 10-5 

2.43 X 10-5 

7r/3 

Case 3 

0.040 

0.320 

1.0 X 10-5 

—7r/4 

0.321 

0.0244 

0.639 

7.54 X 10-5 

2.43 X 10-5 

7r/3 


TABLE I: Input values of the parameters for the three scenarios considered for the four-neutrino hypothesis. Values of 
013, and 023 are chosen to be consistent with the best-fit values of |Ue 2 p, |Ue 3 |^, and |U;j 3 p, given choices of 0i4 and 024. Here, 
rjs = ri 2 — ri 3 . Note that Ami 4 is explicitly assumed to be positive, i.e., ml > mf. 

In all three Cases, we assume that the neutrino mass hierarchy for the mostly active states is normal, i.e. Am |3 = 


t Note, however, that Case 2 is in slight disagreement with existing bounds from Daya Bay and MINOS, see Fig, 
§ A recent analysis of solar and reactor data constrain sin^ 0i4 < 0.04 m- This bound is not depicted in Fig. 2 










































+2.43 X 10“^ eV^, and in all Cases we assume r]i = tt/S and rjs = — 7r/4, typical of scenarios where CP-invariance 
violating effects are large.^ For completeness, we also assume, in all Cases, sin^ (^34 = 0. Gaussian priors are adopted, 
mostly from solar neutrino data and data from KamLAND, on the solar parameters, |Pe 2 p = 0.301 ± 0.015 and 
Amj 2 = (7.54 ± 0.24) x 10“^ eV^ [52]. Without these priors, DUNE is mostly insensitive to either Am 22 or cj)i 2 . We 
make use of the Markov Chain Monte Carlo package emcee |3S] , which estimates a probability distribution for each 
fitting parameter. Figs. 11 1^ and 1^ in Appendix B, depict sensitivity contours at 68.3%, 95%, and 99% CL and 
one-dimensional distributions for the ten parameters, for Cases 1, 2, and 3, respectively. Input values from Table|l] 
are shown as stars in the two-dimensional plots. Given the fact that the amount of information in Figs. [min and 
13 is somewhat overwhelming, in order to guide the following discussions, Figs.|^[^ and [^depict sensitivity contours 
at 68.3%, 95%, and 99% CL for a subset of the parameters of interest, for Cases 1, 2, and 3, respectively. 

Fig. [^depicts the fit results for a subset of the parameters (sin^(/) 24 , sin^(/)i 4 , r]i, rj^ and assuming Case 1. 

Here, the values of sin^ (j )24 and Amf^^ can be excluded from 0 at the 99% CL, while the value of sin^ ^14 is consistent 
with 0 at 68.3% and rjs cannot be constrained at the 95% CL. Nonetheless, the CP-odd phase 771 , which can be 
more or less trivially associated with the CP-odd phase Scp in the three-neutrino scenario, is constrained to be 
nonzero at the 99% CL. As expected, there is very little sensitivity to Am^ 4 , except for establishing that it is large 
(Amf 4 > 7.9 X 10“^ eV^ at the 99% CL). 


10^1 


-e- 


10 


,-2 


10 " 


10 

10 


> 

(D 


<1 10 " 


10 


-2 

10 



sin^ (^24 


Vi 


FIG. 4: Expected sensitivity contours at 68.3% (blue), 95% (orange), and 99% (red) CL at DUNE with six years of data collection 
(3y + 3y F), a 34 kiloton detector, and a 1.2 MW beam given the existence of a fourth neutrino with parameters from Case 1 

Tableln Results from solar neutrino experiments are included here as Gaussian priors for the values of |Ue 2 p = 0.301 ±0.015 
and AmC = 7.54 ± 0.24 x 10"® eV^ 


^ We explored several other sets of input values for rji and rjg. This particular choice leads to generically large effects without extraordinary 
cancellations, enhancements, or ambiguities. 
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Fig. [^depicts the fit results for a subset of the parameters (sin^(^ 24 , sin^(^i 4 , r]i, and r]s), assuming Case 2. Here, 
the values of sin^ (j)i 4 , sin^ (^ 24 , (cf. Fig. [^, and rjs are observed at at least the 95% CL, i.e., the fit establishes 

that none of the new physics parameters vanish. In particular, the values of sin^ (/)i 4 , sin^ ^ 24 j and Amf^ are excluded 
from zero at the 99% CL. In this case, there is enough sensitivity to the two independent CP-odd phases to establish 
that not only there are new neutrino degrees of freedom but that there is more than one new CP-invariance violating 
parameter in the theory. In summary, one can establish that there is new physics beyond the standard paradigm, and 
that the new physics is CP-invariance violating. 



3 years z/ + z/ 

- (7.54 ± 0.24) x 10“^ eV^ 
AmK = 2.43 X 10"^ eV^ 

Amf4 = 10“^ eV^ 
sin^ <l)i2 = 0.315 
sin^ = 0.024 
sin^ <^23 = 0.456 
sin^ <^14 = 0.022 
sin^ <p2i = 0.030 
Vi ^ 

'Hs ~ 

=0.301 ± 0.015 


FIG. 5: Expected sensitivity contours at 68.3% (blue), 95% (orange), and 99% (red) CL at DUNE with six years of data collection 
(3y -I- 3y F), a 34 kiloton detector, and a 1.2 MW beam given the existence of a fourth neutrino with parameters from Case 2 

in Table^ Results from solar neutrino experiments are included here as Gaussian priors for the values of |Pe 2 |^ = 0.301 ±0.015 
and Am^ = 7.54 ± 0.24 x 10"® eV^ 1^ . 


Fig.|§ depicts the fit results for a subset of the parameters (sin^ 4>24, sin^ (j)i 4 , and Amf^), assuming Case 3. The 
results here are somewhat similar to (but less constraining than) those from Case 1. The measurement of Amf^ is 
consistent with 0 at 68.3% CL. but, as expected, the data reveal that it is small {Aml^ < 1.6 x 10“^ eV^ at the 99% 
CL). The new CP-odd phase cannot be measured significantly (cf. Fig. 13). On the other hand, one can exclude the 
hypothesis that the “standard” CP-odd phase rji is zero, but the sensitivity is worse than what one can achieve if the 
data were consistent with the three-flavor scenario. 

If there is a fourth neutrino mass-eigenstate, the parameters of the fourth neutrino may significantly affect DUNE’s 
ability to measure the mixing angles naively associated with three-neutrino oscillation. For example, as shown 
in Appendix A and in Ref. [T], the expected measurement precision for 9i^ assuming a three-neutrino scenario is 
S0i2/9i3 ~ 3%. In Cases 1 and 3, this precision is much worse, S(j)i 3 /(l)i 3 — 10%. In Case 2, however, the precision 
with which (j)i 3 can be measured is S4ii3/(j)i3 ~ 4%, i.e., similar to the precision obtained in the three-neutrino scenario. 
This happens because, in Case 2, one can mostly disentangle effects due to the different oscillation frequencies. 


C. Testing the Three-Massive-Neutrinos Paradigm 


In Sec. |III A[ we simulated data assuming a three-neutrino scenario and, by analyzing it assuming the four-neutrino 
hypothesis, were able to constrain the values of the new mixing parameters. In Sec. |IIIBl we simulated data assuming 
different four-neutrino scenarios and, by analyzing it assuming the four-neutrino hypothesis, were able to constrain 
or measure, sometimes quite precisely, the new mixing parameters. Here we address a different question; if we were 
to simulate data consistent with a four-neutrino scenario, would we be able to tell that there are more than three 
neutrinos? More concretely, would the analysis of the data assuming the three-massive-neutrinos paradigm reveal 

that the paradigm is incorrect? _ 

To address this question, we fit the expected yields from Case 2, introduced in Sec. HIB (see Table|T| assuming the 
three-neutrino hypothesis. We obtain best-fit values of 6 * 12 , ^ 13 , ^ 23 : Aml 2 , Am\^, and 6 cp, along with associated 
uncertainties. The precision with which the parameters can be measured is comparable to what would be expected of 
DUNE if the data were consistent with the three-neutrino hypothesis. We also find, however, that the overall quality 
of the fit is poor: degrees of freedom (dof) ~ 180/114, or a discrepancy of roughly 4cr. Hence, the three-neutrino 
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sin^ 014 


3 years u-l-u 

Am ^2 = (7.54 ± 0.24) x eV^ 
Ami\ = 2.43 X 10“^ eV^ 

Ami^4 = lO-'* eV^ 
sin^ 042 = 0.321 
sin^ 043 = 0.024 
sin^ 023 = 0.639 
sin^ 044 = 0.040 
sin^ 024 = 0.326 

m = V3 
% = 

|;7,2|^ =0.301 ± 0.015 



sin^ 024 


FIG. 6: Expected sensitivity contours at 68.3% (blue), 95% (orange), and 99% (red) CL at DUNE with six years of data collection 
(3y + 3y F), a 34 kiloton detector, and a 1.2 MW beam given the existence of a fourth neutrino with parameters from Case 3 

in Tablellj Results from solar neutrino experiments are included here as Gaussian priors for the values of |f/e 2 p = 0.301 ±0.015 
and Am^ = 7.54 ± 0.24 x 10"® eV^ 


hypothesis cannot mimic the additional oscillations associated with Atoi 4 = 10 “^ eV^, for any set of values of the 
three-neutrino parameters.** 

Once a bad goodness-of-fit is established, it becomes crucial to identify in which way the three-neutrino hypothesis 
fails. This can be done in a variety of ways. Here, for illustrative purposes, we try to diagnose the poor goodness-of-fit 
by splitting the data set into two subsets: the appearance data and the disappearance data, and analyze both subsets 
separately (combining neutrino and antineutrino data in each case). In both sub-channels, the extraction of 612 and 
Am 42 is mostly driven by the priors from solar neutrino data, while the disappearance data are mostly insensitive to 
the CP-odd parameter Scp- For these two reasons, it is most illuminating to examine the measurements obtained from 
these two fits in the sin^ ^13 - sin^ 623 plane, depicted in Fig. Fig. [^reveals that the appearance and disappearance 
channels favor different values of sin^ 9 i 3 and sin^ 623 , with no overlap of the preferred regions at the 68.3% CL. For 
the appearance channels, the fit has Xmin/ ^of ~ 78/54 (roughly 2a), and for the disappearance channels, Xmin/ 

~ 91/54 (roughly 3a). The overall four-sigma (very significant) discrepancy, therefore, is, in some sense, the product 
of a mediocre fit in the appearance channel, a poor fit in the disappearance channel, and the fact that the two subsets 
of data point to different regions of the parameter space. 

The shapes observed in Fig.j^are easy to understand. In a three-neutrino scheme, the disappearance probability 


While we concentrate on Case 2 here, we obtain poor fits also by assuming data consistent with Cases 1 and 3, where the new oscillation 
frequency cannot be explicitly observed. The discrepancy is most significant for Case 2, however. 
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FIG. 7: Expected sensitivity contours at 68.3%, 95%, and 99% for neutrino and antineutrino appearance channels (blue, orange, 
red), and neutrino and antineutrino disappearance channels (green, teal, blue) in the sin^ - sin^ 023 plane, assuming the 
data are consistent with Case 2 and analyzing it assuming the three-massive neutrinos paradigm. 


depends mostly on the |I 7 e 3 p(l— The best fit values translates into the relations sin^ Ois ~ 0.43(l-|-sin^ ^ 13 ) or 

sin^ 6*13 0.57(l-|-sin^ 6*13), which explains the approximately linear shapes in Fig.j^ The appearance channels, on the 

other hand, are mostly sensitive to the product sin^ 6*13 sin^ 623 , which explains the hyperbolic shape in Fig. Hence, 
in order to fit the four-neutrino data, the two different data sets wander towards different best-fit values for sin^ 0 i 3 
and sin^ 623 as they strive to maintain the different combinations of these parameters constant. In Appendix A, we 
repeat this analysis, this time simulating data consistent with the three-neutrino scenario. The results are depicted 
in Fig. The shapes obtained from the two subsets are similar to those in Fig. but in this case the two analyses 
point to the same values of sin^ ^13 and sin^ 623 . 

In summary, not only is the goodness-of-fit poor, it is also possible to ascertain that different measurements of 
the mixing parameters are inconsistent with one another if one assumes that the three-massive-neutrinos paradigm 
is correct. There are several other “inconsistency checks” one would perform in order to reveal that new physics is 
affecting the long-baseline oscillations, including comparing data obtained with the neutrino beam and the antineutrino 
beam, comparing DUNE data with those from HyperKamiokande (same LjE^ values, but different neutrino energies 
and baselines), comparing DUNE data with data from “ 0 i 3 ” reactor neutrino experiments |36H38j . medium baseline 
reactor experiments [39l[40], atmospheric neutrino experiments (for example, PINGU |4T]), etc. 


IV. SUMMARY AND CONCLUSIONS 

Very ambitious next-generation long-baseline neutrino oscillation experiments are currently under serious consid¬ 
eration, especially the “superbeam” experiments Fermilab to DUNE in the United States and J-PARC to Hyper¬ 
Kamiokande in Japan. Among the goals of these projects are searching for CP-invariance violation in the lepton 
sector and testing the limits of the three-massive-neutrinos paradigm. Here, we addressed the capabilities of the 
DUNE experiment to discover a fourth neutrino mass-eigenstate or, instead, constrain its existence, either falsifying 
or strengthening the three-massive-neutrinos paradigm. While several different new phenomena could manifest them¬ 
selves at long-baseline neutrino experiments, we chose one new neutrino mass-eigenstate for a few reasons. First, 
oscillation effects due to a new light neutrino mass-eigenstate are easy to parameterize, and very familiar. Second, 
light sterile neutrinos are a very natural and benign extension of the Standard Model and could, for example, be a side 
effect of the mechanism responsible for the nonzero neutrino masses. Finally, the so-called short-baseline anomalies 
may be interpreted as evidence for new neutrino degrees of freedom so it is possible, even though the evidence is not 
very robust, that new neutrino states have already been found. 
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Assuming coherent oscillations, we discuss the oscillation probabilities involving a fourth neutrino for a wide range 
of values for the new mass-squared difference, including |Amj 4 | 3> |Am 43 |, when the new oscillation length is very 
short and expected to lead to averaged-out effects at DUNE, \Aml^\ lAmfgj, when the new oscillation length is 
too long to be observed at DUNE, or |Amf 4 | ^ jAm^gl, when DUNE is sensitive to the new oscillation frequency. We 
highlight the fact that, in all three cases, the values of the active elements of the fourth column of the mixing matrix, 
Uq, 4 , Of = e, /i, T, have a nontrivial impact on the experiment as long as these are large enough. We also discuss the 
extra sources of CP-invariance violation that arise from phases in the new mixing matrix elements. Given access to 
the r'e-appearance and t'^-disappearance channels, we find that DUNE is, in principle, sensitive to two of the three 
CP-odd phases in the mixing matrix. 

We simulate data in the DUNE experiment assuming a 34 kt detector, a 1.2 MW proton beam, and 3 years each 
of neutrino and antineutrino data collection, exploring different scenarios. If the data are consistent with three- 
neutrinos (i.e., there are no accessible new light neutrinos) we find that DUNE is less sensitive than, for instance, 
the Daya Bay experiment when it comes to constraining |Pe 4 p if I ^”^ 14 ! ^ while DUNE can outperform 

current long-baseline experiments when it comes to constraining |P/j 4 p if |Ami 4 | > lAmfgj. On the other hand, 
if |Atoj 4 | jAm^gl, DUNE outperforms all current experiments when it comes to constraining new, light neutrino 
mass-eigenstates thanks, in part, to the broad-band-beam nature of the experiment and the fact that it measures 
both 14^ disappearance and Vg appearance. 

If the data are consistent with the existence of a fourth neutrino, DUNE has the capability to measure the new 
mixing parameters. This capability, however, depends strongly on the values of the parameters associated with the 
fourth neutrino, particularly Aml^. We find that there are circumstances under which DUNE can not only discover 
new physics but also establish that there are new sources of CP-invariance violation. We emphasize that, if there is 
a new neutrino mass-eigenstate, the Vg and data at DUNE can only explore a subset of the existing parameter 
space. One of the new mixing angles, and one of the two new sources of CP-invariance violation can only be accessed 
if one could also study Vr appearance (or construct a j^T--beam, a much more challenging proposition). 

Improved sensitivity is expected, of course, if backgrounds at DUNE turned out to be smaller than anticipated, if 
better energy resolution can be achieved, or if one were to optimize the beam energy profile. For example, we find 
that a 50% background reduction yields C(10%) improved sensitivity to a new neutrino state (e.g., C(10%) stronger 
constraints in the 4 |Pe 4 p|P^ 4 p - Am\/^ plane, cf. Fig. Access to higher neutrino energies, on the other hand, 
would allow DUNE more sensitivity to higher values of Ato^ 4 . 

We also briefly addressed whether DUNE data could reveal the existence of physics beyond the three-massive- 
neutrinos paradigm if the data were consistent with the existence of a fourth neutrino. We find that DUNE data are 
precise enough to reveal that a three-neutrino fit to data consistent with four neutrinos is poor (assuming the new 
mixing parameters are accessible). We also show that, in this scenario, fits to disjoint subsets of DUNE data point to 
different regions of the three-neutrino parameter space, another sign of new physics beyond three active neutrinos with 
nonzero mass. In order to properly diagnose that (a) there is physics beyond the three-massive-neutrinos paradigm, 
and (b) determine the nature of the new physics, it is very likely that one will need more/better data. We hope to 
return to these very important issues in a future study of long-baseline neutrino oscillations and new phenomena. 
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Appendix A: Three-neutrino Fits to Three-neutrino Data 

In this Appendix we present the results of simulating and analyzing DUNE data consistent with the three-neutrino 
hypothesis. This is done, in part, to validate our simulation and analysis tools, and in order to facilitate comparisons 
between the three-neutrino and the four-neutrino scenarios. We also comment on the ability of the DUNE experiment 
to constrain the solar parameters Amf 2 and 612 . 

Fig. 1^ depicts the expected sensitivity contours at 68.3% (blue), 95% (orange), and 99% (red) CL as measured 
by DUNE with six years of data collection (3 years with the neutrino beam, three years with the antineutrino one), 
a 34 kiloton detector, and a 1.2 MW beam, assuming the data are consistent with a three-neutrino scenario. The 
figure also depicts one-dimensional A^ plots for each parameter, with the 68.3% (blue), 95% (orange), and 99% 
(red) CL highlighted. Quoted measurement bounds are for 68.3% CL. The input values of the mixing angles are 
sin^ 012 = 0.308, sin^ 0 i 3 = 0.0235, and sin^ 023 = 0 .437. Results from solar neutrino experiments and KamLAND are 
included here as Gaussian priors for the values of |Ue 2 p = 0.301 ± 0.015 and Am \2 = (7.54 ± 0.24) x 10“^ eV^. This 
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scenario assumes a normal mass hierarchy, i.e., = +2.43 x 10“^ eV^, and that the CP-odd phase is ^cp = tt/S. 

Distributions are sampled using a Markov Chain Monte Carlo (MCMC) method [35]. We find that the measurement 
precisions of the mixing angles and mass-squared differences are safely comparable to the projected results in Ref. [T|. 



sin ^13 — 0-023^Qom 
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sin^ 0^2 = 0.308 
sin^ 0^3 = 0.023 
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Scp = 1-047 
\UJ'^ =0.301 ± 0.015 




FIG. 8 : Expected sensitivity contours at 68.3% (blue), 95% (orange), and 99% (red) CL as measured by DUNE with six 
years of data collection (3y p + 3y u), a 34 kiloton detector, and a 1.2 MW beam given a three-neutrino scenario. On the 
far-right, one-dimensional plots for each parameter display the 68.3% (blue), 95% (orange), and 99% (red) CL bounds. 
Quoted measurement bounds are for 68.3% CL. Mixing angles here are sin^ 612 = 0.308, sin^ ^13 = 0.0235, and sin^ 62 i = 0.437. 
Results from solar neutrino experiments are included here as Gaussian priors for the values of \Ue 2 \^ = 0.301 ± 0.015 and 
Ami 2 = (7.54 ± 0.24) x 10“® eV^. This scenario assumes the normal hierarchy, i.e. Arrii^ = +2.43 x 10“^ eV^, and that 
Scp = p/3. 


Throughout, with one exception, our analyses include only data to be collected by the DUNE experiment. We do 
not include, for example, data from “^la” reactor experiments, [36H38j nor do we include existing or simulated data 
from the long-baseline experiments currently in operation [42H45j . The main reason for this is that we anticipate 













































































































14 


DUNE data to provide the most significant information when it comes to measurements of sin^ 613 , sin^ 023, ^^13 
(including the sign), and Scp- We do, however, include results from solar data and from the KamLAND experiment 
when it comes to constraining the solar parameters Amf 2 and 0 i 2 - The reason for this is that DUNE’s ability to, in 
isolation, determine the solar parameters is rather limited. To illustrate this fact, we perform a DUNE-only ht to the 
DUNE data. The results for a subset of the three-neutrino oscillation parameters (Am 22 , tan^ 0 i 2 , Scp) are depicted 
in Fig. 
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FIG. 9: Expected sensitivity contours at 68.3% (blue), 95% (orange), and 99% (red) as measured by DUNE with six years of 
data collection (3y p + 3y V), a 34 kiloton detector, and a 1.2 MW beam assuming no fourth neutrino, and without external 
information from solar neutrino experiments included as Gaussian priors. 


While DUNE is not very sensitive to tan^ 0 i 2 or Ami 2 (note the logarithmic scales), it can exclude nonzero values 
for each parameter, and is still able to observe CP-invariance violation and measure Scp even if external information 
on the solar parameters is not included in the data analysis. The uncertainty on Scp, as expected, is significantly 
larger (cf. Fig. [^. 

Finally, we repeat the analysis discussed in Sec. [Ill C| where appearance and disappearance channels are analyzed 
independ ently, t his time assuming the data are consistent with the three-neutrino scenario. Fig. |10| depicts the results. 
As in Sec. Ill C we see that the appearance channel is sensitive to the product sin^ 0 i 3 s in^ 023 while the disappearance 
channel is mostly sensitive to |17^3p(l — |t/^ 3 p). Unlike the scenario discussed in Sec. IIIC (see Fig. [7|, here the fits 
to the different data sets are in agreement. 
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FIG. 10: Expected sensitivity contours at 68.3%, 95%, and 99% for neutrino and antineutrino appearance channels (blue, 
orange, red) vs. neutrino and antineutrino disappearance channels (green, teal, blue) in the sin^ dn - sin^ 623 plane assuming 
a three neutrino hypothesis with parameters from Ref. [^, indicated with a star in the figure. 


Appendix B: Four-neutrino Fits to Four-neutrino Data 


Here we display the full results we obtain when analyzing the different four-neutrino scenarios (see Table M assuming 
the four-neutrino hypothesis. Figs. 11 KHl and [I^ depict the expected sensitivity contours at 68.3% (blue), 95% 
(orange), and 99% (red) CL at DUNE with six years of data collection (3y i' -b 3y F), a 34 kiloton detector, and 
a 1.2 MW beam, given the existence of a fourth neutrino with parameters from Case 1, Case 2, and Case 3 in 
Table |lj respectively. Results from solar neutrino experiments are included here as Gaussian priors for the values of 
|t^e 2 p = 0.301 ± 0.015 and Ami 2 = 7.54 ± 0.24 x 10“® eV^ |22]. Distributions are sampled using a Markov Chain 
Monte Carlo (MCMC) method [3S]. These results are discussed in Sec. IIIB 
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FIG. 11: Expected sensitivity contours at 68.3% (blue), 95% (orange), and 99% (red) CL at DUNE with six years of data 
collection (3y z/ + 3y F), a 34 kiloton detector, and a 1.2 MW beam given the existence of a fourth neutrino with parameters 
from Case 1 in Table On the far right, one-dimensional plots for each parameter display 68.3% (blue), 95% (orange), and 
99% (red) CL bounds. Quoted measurement bounds are for 68.3% CL. Results from solar neutrino experiments are included 
here as Gaussian priors for the values of \Ue 2 \^ = 0.301 ± 0.015 and Amf 2 = 7.54 ± 0.24 x 10~® eV^ |22| . 
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FIG. 12: Expected sensitivity contours at 68.3% (blue), 95% (orange), and 99% (red) CL at DUNE with six years of data 
collection (3y i/ + 3y F), a 34 kiloton detector, and a 1.2 MW beam given the existence of a fourth neutrino with parameters 
from Case 2 in Table On the far right, one-dimensional Ax^ plots for each parameter display 68.3% (blue), 95% (orange), and 
99% (red) CL bounds. Quoted measurement bounds are for 68.3% CL. Results from solar neutrino experiments are included 
here as Gaussian priors for the values of \Ue 2 \'^ = 0.301 ± 0.015 and Am \2 = 7.54 ± 0.24 x 10~® eV^ |22| . 
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sin^ <^12- 0.31+0*1 


3 years iy-\-iy 

Am^2 = ( 7-540 ± 0 . 240 ) X 10 “'^ eV^ 
Ami ^3 = 2.430 X 10 “^ eV^ 

Ami\ = lO-'^ eV^ 
sin^ (pi2 ~ 0.321 
sin^ <)ii 3 = 0.024 
sin^ (l>23 = 0.639 
sin^ 4>u = 0.040 
sm" 024 = 0.326 
7?i = 7r/3 
Vs ~ ~'^/4 
\U,2'? = 0.301 ± 0.015 



+0.020 win-S „-ir2 


sin'^ ^12 


sin <5i>23 


sin^ ^14 


Ami2[10 ^eV^]Ami3[10 ^ eV^ ] Am44[eV^] 


FIG. 13: Expected sensitivity contours at 68.3% (blue), 95% (orange), and 99% (red) CL at DUNE with six years of data 
collection (3y !/ + 3y V), a 34 kiloton detector, and a 1.2 MW beam given the existence of a fourth neutrino with parameters 
from Case 3 in Table On the far right, one-dimensional plots for each parameter display 68.3% (blue), 95% (orange), and 
99% (red) CL bounds. Quoted measurement bounds are for 68.3% CL. Results from solar neutrino experiments are included 
here as Gaussian priors for the values of \Ue 2 \^ = 0.301 ± 0.015 and Amf 2 = 7.54 ± 0.24 x 10~® eV^ |22| . 


























































































































































































































































































































